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Abstract We develop a unified approach to integrating the Whitham modulation equa- 
tions. Our approach is based on the formulation of the initial value problem for the zero 
dispersion KdV as the steepest descent for the scalar Riemann-Hilbert problem and on 
the method of generating differentials for the KdV- Whitham hierarchy 0. By assuming the 
hyperbolicity of the zero-dispersion limit for the KdV with general initial data, we bypass 
the inverse scattering transform and produce the symmetric system of algebraic equations 
describing motion of the modulation parameters plus the system of inequalities determining 
the number the oscillating phases at any fixed point on the x,t - plane. The resulting system 
effectively solves the zero dispersion KdV with an arbitrary initial data. 

1 Introduction 

The initial value problem for the Korteweg - de Vries equation 

ut - Quur, + e^Uxxx = 0, u{x, 0) = uo{x) (1) 

in the zero dispersion limit e — > has been the subject of intense studies during more than 
25 years. The physical interest of this limit is that it allows to model the phenomenon 
of dispersive shocking in dissipationless dispersive media. In contrast to the usual dissi- 
pative hydrodynamics, the regularization of a shock occurs here through the generation of 
small-scale nonlinear oscillations. Gurevich and Pitaevskii in [Q] proposed to describe these 
oscillations with the aid of the one-phase Whitham modulation equations Multiphase 
analog of the Whitham equations was derived by Flaschka, Forest and McLaughlin ||^ . Later, 
Lax and Levermore 0] rigorously showed that the multiphased averaged equations appear in 
the zero dispersion limit of the initial value problem for the KdV equation with asymptot- 
ically reflectionless initial data. For the reflecting potentials the zero dispersion theory was 
constructed by Venakides who also identified the parameters describing the weak limit of 
the solution in the Lax-Levermore approach with the Riemann invariants of the modulation 
equations in p (see [|], 



The modulation system itself, however, has to be integrated to reveal dependence of the 
Riemann invariants (modulation parameters) on x and t. For this, the observation of Tsarev 



(see also [11 ) who generalized the classical hodograph method to a multidimensional (in 
the space of dependent variables) case was crucial. The result of applying this generalized 
hodograph transform to the averaged KdV is the overdetermined consistent system of linear 



PDE's. Tsarev's results were put into the algebro-geometrical setting by Krichever |12|, [13 



Each solution to the Tsarev system gives rise to some local solution of the nonlinear 
modulation system. These solutions generically exist only within definite regions of the x, t 
- plane. To obtain the global solution one should supplement the constructed local solutions 
with the information about the number of nonlinear phases in each region and provide 
the smooth matching of solutions with different N on the phase transition boundaries. The 
simplest, yet very important, class of problems with < 1 has been investigated in works of 



Tian [14[, [|T^, Gurevich, Krylov, El and their collaborators [1G[, [17] who obtained a number 
of exact solutions for the initial value problems with monotone and hump-like data having 
the only break point. Tian [jTSf was also able to prove that the solutions to the Whitham 
system in this case do globally belong to iV = 1. 

More general local solutions to the modulation system in the case of arbitrary were 
constructed in a symmetric form by El [Q with the aid of fundamental solution (an analog 
of the Green function) to the Tsarev equations. 

Some global solutions involving the case N = 2 were constructed recently by Grava |T^ 
on the basis of Dubrovin's variational approach to the Whitham equations [p!9| . 

The common feature of all pointed methods for constructing the global solutions to the 
Whitham equations is the necessity to follow all earlier times t < to to obtain the solution at 
to. This difficulty has been bypassed resently by Deift, Venakides and Zhou [0] who obtained 
the algebraic equations for the KdV modulations avoiding the integration of the modulation 
equations themselves. The method of DVZ is based on reformulation of the initial value 
problem for the zero dispersion KdV with decaying solitonless analytic initial data as the 
steepest descent for the scalar Riemann-Hilbert problem for the complex phase function. As 
a result, they obtain not only the local solutions to the modulation equations, but also the 
system of inequalities enabling to determine the number of phases locally at each point of 
the X, t-plane. 

It is clear, however, that the requirements of decaying and pure reflection for the initial 
potential are not essential in the zero dispersion limit. Really, assuming the finite speed of 
propagation (i.e. validity of the Whitham equations) for the zero-dispersion KdV with an 
arbitrary initial data one can see that only the finite part of the initial data contributes to the 
solution at any finite t. In this case, considering the evolution of 'multihump' potential and 
and then tending number of 'humps' to infinity one arrives at the solution for nondecaying 
initial data. Certainly, for the finite e the results obtained in this way are valid within finite 
time interval until the semiclassical (Whitham) spectrum begins to compete with the fine 
spectrum of the potential. The assumption of hyperbolicity in the limit studied allows also 
to omit the requirement of analyticity for the initial data. 

In this work, we combine methods of Deift, Venakides, and Zhou ^ and El [§ to produce 
the global solution to the KdV modulation equations with an arbitrary initial data. Namely, 
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by assuming validity of the modulation equations in the case of the zero dispersion KdV with 
an arbitrary initial data we reformulate the initial value problem for the Whitham system 
as the Riemann-Hilbert problem and produce the symmetric system of algebraic equations 
supplemented by the system of inequalities determining both the motion of the Riemann 
invariants and the change of the genus of the Riemann surface (number of oscillating phases). 
We represent the local part of the obtained general solution in a potential form with the 

potential. Recently the 



generalized functional of the Peierls-Frohlich type I^O] , I^T] , 
functionals of this type in more particular form were shown to give rise by minimization to 
the global solutions for the Whitham equations with monotone initial data | |19| ] . 

We also make identification of the part of the results obtained in [Q] with the results ob- 
tained earlier in [0 , and [|l^] which unifies different approaches developed independently 
in this area. 



2 Summary of the Riemann-Hilbert Steepest Descent 
Method for the Zero Dispersion KdV with Decaying 
Solitonless Initial Data 

It is well known that the inverse scattering transform can be reformulated as a matrix 
Riemann-Hilbert (RH) problem (see [|^, [|4|)- Deift, Venakides, and Zhou showed by 



considering the quasiclassical asymptotics for the initial value problem (|I|) in the case of 
one-hump solitonless initial data < Uq{x) < 1 that this problem can be asymptotically (as 
e — *• 0) reduced to the scalar RH problem for the complex phase function g{X). We present 
here the resulting formulas (for details see [0, p5||). 



Let the interval (0, 1) of the real axis on the A - plane is partitioned into a finite set of 
intervals {Gk{k = 0, 1, . . . , A^) : (0, n), (rs, rg), . . . , (ra^v, Tstv+i)}, {fi„(n = 1, . . . , + 1) : 
{rur2), (r3,r4), . . . , {r2N+i, !)}• 

We introduce the combination 



a 



{X) = 4tX^/^ + xX'/\ (2) 



which will be important hereafter, and the phases of the reflection (from the right) and the 
transmission coefficients for the scattering on the potential uo{x). 

oo 

P+(A) = X+Al/2 ^ J [^1/2 _ _ ^^(3;'))l/2]^^/ ^ (3) 
x+ 
x+ 

r(A) = [ iuoix') - Xy/^dx' . (4) 



where = a;^(A) are the roots of the equation uo{x) = X. 

Then the RH problem defining the complex phase function g{X), which plays the central role 
in the future, takes the form 
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On the intervals {Gk, {k = 0, . . . , N) : (0, ri), (r2, r^), . . . , {r2N, r2N+i)} the following rela- 
tions hold: 

^ = p;-a' (5) 

-.<^<0. (6) 

It follows from the equality (||) that on each interval Gk- g+ + g- — 2p+ + 2a = —Qk, 
where Qk is some constant of integration; without loss of generality one can put Qq = 0. 
On the intervals {Bk, {k = 1, . . . , N) : (ri, r2), {r^, ^4), . . . , {r2N-i,r2N)} we have 

(7) 

<P+-a . (8) 

The latter ineguality is due to the further reduction of the RH problem with the aid of the 
steepest descent method |l|]. 

On the remaining interval {-Bw+i : {r2N+i, 1)} there exist two possibilities: 

A 9+ -9- 9'+ + 9'- ^ , I /ox 

case A: — - — = -r , < p^[\) - a . (9) 

which coincides with the conditions ([7p, (|p 
or 

caseB: g^-g_ = Q^ ^+±^ > p'^(A) - . (10) 

Also we have outside the interval (0, 1) 

g+ + g- = ^ if A<0, (11) 

g+-g-=Q if A>1. (12) 
In all above formulas we denoted 

g±{\) = \\mg{\±i6) 
The additional requirement imposed on the function g{X) is 

The functions {^J~Xg' {\))± are continuous for real A , (13) 
Also it follows from |]l| that the following asymptotic holds 

g{\)=g^/X^I^ + 0{l/X) as A ^00. (14) 
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One can see that just formulated RH problem contains not only equations defining the 
function g{X) on the intervals Bj. and but also the inequalities which define the overall 
number 2N of the intervals (if N has been chosen wrongly then at least one of the inequalities 
will be violated). 

We observe that for each fixed x and t, g' also satisfies a scalar RH problem on the real 
axis. Indeed, g'j^ + g'_ = Q when A < 0, and g'_^ — g'_ = when A > 1. When < A < 1, 
the equalities (^, (|^), (^ and (|10|) specify g'^ + gL — 2p' + 2a' = on each interval Gk, 
while g'j^ — g'_ either equals or equals — 2zr' outside these intervals. We only consider 
g'^ — gL = — 2zr' when A lies between any two of the G^s while on the remaining interval 
{r2N+i, 1)5 we examine both possibilities i.e. g'^ — gL = —2iT' (case A) or g'_^_ — g'_ = (case 
B). 

Necessarily g'{X) has the following form. 



2p'(/i) - 2a'{fx) dfi f -2zr'(/i) dfi 



(15) 



where 



and E = 1 in the case A and E = r2N+i in the case B. Here, \/-R2iv+i(A) is positive for 



2Af+l 

i?27v+i(A)= n (^-^j)' (16) 



A > 'r2Ar+i. Also, y RtN+iW denotes the boundary value from above. 

A sufficient number of conditions to determine the endpoints of the intervals Gj can now 
be written down. Indeed, the condition g{X) = 0(A~^/^) for large A implies g'{\) = 0(A^^/^), 
which leads to the following moment conditions, 

- X^dX + I r'"'^^^ XHX = , (17) 



k = 0,...,N. 

A second set of conditions is obtained by integrating g' around each Gj and using (0): 
In case A we obtain 



G, 



(g'iX) - g'_iX))dX = -2z(r(r2,+i) - r(r2,)), j = 1, . . . AT . 



In case B, r(r2Ar) and T{r2N+i) niust be replaced by zero. 

Conditions ([T7|) and ([T8| ) represent a system of {N + 1) + = 2A^ + 1 independent 
equations for the 2A^ + 1 unknowns (the branch points ri, . . . , r2N+i of the Riemann surface 
(p!6|) ). Conversely suppose that for given x, t and some A^, the quantities ri, . . . , r2Ar+i satisfy 
conditions { \L7\ ) and (p!8D, giving rise to an explicit expression (|1^) for g', and hence for g by 
integration. Suppose further that the function g so constructed, also satisfies the inequahties 
(^ - (0). Then g is the desired solution of the scalar RH problem. 
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Deift, Venakides and Zhou show that the parameters rj represent the "semiclassical 
spectrum" of the problem, i.e. they are the branch points of the spectral Riemann surface 
defining the local finite-gap solution of the KdV. In Sec. 4 we shall directly identify the 



moment conditions ([T7|) , ([T8|) with the local solution of the KdV-Whitham system. 

Now we present a number of important relations which will be very useful for the iden- 
tification. 



We observe from the relations (Bf) - (IT^), that g{X) also satisfies a RH problem. Solving 



this RH problem in exactly the same way as the problem for g' we obtain similarly to (15) 



j=0,...,N 



2p^{fi) — 2a{jj,) — ^Ij djj, r — 2ir(/i) djj, 

(19) 

The expression for the constant of integration Qj is 

rij = -2x <f -8t<f X^'^i)j + 4 / p+^ljj - 4 / iTi)j . (20) 

JuGk J{0,E)\uGk 

= xn^i + tnk2 + ^k3 J = i,...,N , no = o. 

We recall that 

E = 1 (case A) or E = r2N+i (case B) (21) 
The basis of holomorphic differentials ipj is given by 

N-l 

= E CfcJ /„ , , (22) 

k=o ^R2N+i[r,X) 

i^j = 5jk, k,j = l,...,N. (23) 

The contours ak {k = 0,1, N) surround the intervals (r2, r^), . . . , (r2j, T2j+i), . . . , {r2N, r2N+i] 
clockwise; 

The identities hold for Qji and Qj2'- 

rij^ = -Res^\^'^i)j , nj2 = -4i?eSooA^/Vi , (24) 
dx^j = ilji , dtVLj = nj2 . (25) 



One can observe from (24) that in normalization (^) accepted, ilji and Q,j2 can be 
identified with the wave number kj and the frequency ujj correspondingly where j is 
the number of the phase , j = 1, . . . ,N Really (taking into account the change of the 
normalization in comparison with 0), we have the expansion as A — >■ cx) 

^^•=A^(^^ + ^^i + •••); (26) 
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where 

2N+1 

kj = CN-l,j, UJj = 2CN-I,j rm + 4:CN-2,j ■ (27) 

m=l 

Taking the mixed derivatives of (^) we arrive at equations expressing the wave number 
conservation laws, for each of phases. 

dtkj = d,uj ,j = l,...,N. (28) 



3 Equations for the Branch Points in a Symmetric 
Form. Identification with Local Solutions to the 
Whitham Equations 

Now we obtain equations for the branch points rj of the Riemann surface, equivalent to the 
moment conditions ([T7| ) , ([T8| ) , in a symmetric form. 
Calculating d^^j directly from (|20|) we obtain 

d,nj = n,i+ -TT^d^r-k, j = l,...,N, (29) 

k=l 



which, together with (^), implies 



2^+^ dVLi 

for any solution rk{x,t). 
Thus, 

djQi = , (30) 

dj = ^, j = l,...,2iV + l, 2 = 1,. ..,iV, 

provided dxTj 7^ , j = 1, . . . , 2N + 1. 

The system (BUf) is the system of A^(A^ + 1) algebraic equations for 2N + 1 variables 



rj{x, t). Due to the uniqueness of the solution for the Riemann- Hilbert problem (^ - (|T^) it 
has to be equivalent to the moment conditions (IT^) , (|l8l) . To show that it is enough to prove 
the consistency of the system (|30D . In other words one has to show that all N closed systems 
for rj{j = 1, 2N + 1) forming the overdetermined system (|30|) are equivalent which is to 
say that (^), in fact, does not depend on the index i. 

With this aim in view we present the following lemma. 
Lemma (Consistency) 

The overdetermined system djQi = , j = 1, . . . , 2N +1, i = 1, . . . , N , where f2j is 
defined by ( PD| ) is consistent and equivalent to the symmetric system of 2N + 1 algebraic 
equations with respect to 2 + 1 variables rj : 
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j{x\^'^ + At\^'^]Kj= j p+{X)Kj~i j r(A)Aj . (31) 

The contours afc, (fc = 0, . . . , A^) as it was mentioned above surround the intervals Gk 
clockwise. The contours /5„, {n = 1, . . . , + 1) surround the intervals i?„ clockwise, and 
the contour Pe surrounds the interval {E, 1), where E = 1 (case A) or = r2Ar+i (case B). 

The differential A^, j = 1, . . . , 2N + 1 is defined by 

N 

~ (A-r,)v/KW ' ^ ^ 

and the coefficients pkj can be found unambigously from the normalization conditions 

'"Aj = 0, m = l,...,N, (33) 



The proof of the lemma can be found in Appendix 1. 

Theorem (Identification) The system of algebraic equations ( PT| ) implicitely defining po- 
sitions of the branch points rj{x, t) of the hyperelliptic Riemann surface in the zero dispersion 
limit of the RH problem for the KdV equation satisfies the A^-phase averaged Whitham-KdV 
system [|], §. 

dtr, = V,ir)d^r^ , j = 1, . . . ,2N + 1 , (34) 

where Vj{r) are computed as certain combinations of complete hyperelliptic integrals [Q. 
Proof 

To identify algebraic system ( pT] ) with the solution of the Whitham equations (Q) we 
make use of the Tsarev result [pT] : 

If Wj{ri, . . . ,r2N+i) is a solution of the linear overdetermined consistent system 

_ d.V, . ^3^^ 



- W, Vi - Vj ' 

then the system of algebraic equations (generalized hodograph transform) 

X + Vjt = Wj (36) 

gives implicitely the smooth solution rj{x, t) to the Whitham system (|5^) provided dxVj ^ 0. 
First, we observe that the following identities hold 
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-TT- f ^^^^^j = 1 for all J , (37) 

A f X'/% = V, . (38) 

The first equality immediately follows from the definition (132|) . The second one requires 
a bit more detailed consideration. We make use of the fact that the system ( [5^ ) implies exis- 
tence of the wave number conservation law (p8[) . Then, introducing the Riemann invariants 
Tj into the equation ( PHD explicitely one easily gets the representation for the characteristic 
speeds of the Whitham system ( p4D (see [0 for the case = 1 and [0] for arbitrary 



N). 

V^ = ^ for any J. (39) 

The expression ( p9D can be interpreted as a generalization of the group velocity notion to 
the case of nonlinear waves @]. It should be noted that the relationships analogous to (39) 
arise in the classical theory of hyperbolic systems as the compatibility condition providing 
the existence of an additional conservation law Using (p^ , (^6]) , ( p2D we arrive at 

TTl Oikj J TTl 

The solution ([3l|) then takes the form (BO) provided 



V, = ±^ f ij^X'/^dX =4: f A^^/^A, = -4pi, . (40) 



p^{X)K-i j r(A)A.}, (41) 



2m 

U«fc U/3„\/3 



where the functions p+(A) , r(A) are supposed to be analytic. 

We have to prove, therefore, that (^) does solve the linear system ( P5| ) which implies 
proving that Aj does solve this system at any A. With this aim in view we consider the 
following combinations occuring in the left-hand part of (|35|) 



where all a^^j are uniquiely defined by the conditions following from the normalization ( p3D 

(9iAj = 0, m=l,...,Ar. (43) 



Another relevant combination is 



A" + E A"-'6i 
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where, again, the coefficients bkj are given by the normahzation 

f{A,-A,) = 0, z^j, (45) 

that imphes bk,i = ak,i- Then 



Ai - Aj pi^i - pij 

Recalhng that pi j = ~l/AVj (see (|40|)) we prove the identification theorem. 
Q.E.D. 

We also present the equivalent form of the solution (|3TD parametrized by the phase of 
the reflection coefficient from the left p_ (cf. (^), 

x~ 

p_(A) = / _ _ uo{x)y^yx - X^'^x- (47) 



It can be shown that the following identity holds (for details see Appendix 2) 



{p_(A) + p+(A)}A, - I f r(A)A, = . (48) 
Then the solution (IsTf) takes the form 



ar 



j {xA^/2 ^ 4tA3/2}Aj- = - ^ p„ {X)Aj + i j T{X)Aj , (49) 

(loo Udfc /3g 

In particular, in the case A we have especially simple representation (we note that Uofc = 
) 

^{a;A^/2 + 4tA=^/2_^p„(A)}Aj = 0, j = 1, . . . , 2Ar + 1 . (50) 

This form of the solution coincides with the one obtained in 0]. 

The differential A^ for the flrst time was introduced by El p| as a generating differential 
of the Whitham hierarchy and can be regarded as a nonlinear analog of the Green function 
for the modulation equations. Really, as Aj depends on a free parameter A and satisfles 
identically the Tsarev equations (P^D, we have proved automatically that it is a fundamental 
solution to the Tsarev system . Expanding (|32|) in powers of 1/A as A oo we obtain 
the homogeneous solutions to Tsarev's equations. Those of them with the odd indices of 
homegeneity n = 3,5, 7. . . are the characteristic speeds of the A^-gap averaged n-th KdV 
in the hierarchy |T^. Therefore, A^ is the generating differential for the averaged KdV 
hierarchy. The effective formulae for the characteristic speeds of the hierarchy in terms of 
the generating differential are (up to a norming constant): 

iy(")(r)dA = -i?eSooA"-3/'A,- , (51) 
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Corresponding solutions rj (x, t) given by the system 

^ {x\'/^ + 4tA^/2 + cA"-2/2}A, = , J = 1,...,2N + 1. (52) 

are self-similar |jl2| 

X 1 

r^-(x,t) =t%-(^), 7= -— ^, n = 3,5,..., (53) 



and solve the initial value problem 



(2n - 1 V 

x=[2«^ p^cw"],=o. (54) 

nl 



Taking = 3 we arrive at the Potemin solution ||2^ (see also |TT[) describing the universal 
regime of the formation of collisionless shock in the vicinity of the break point 0, [|lT|, |1T6 |. 



This solution was also obtained by Wright applying the Lax-Levermore methods [Q. We 
also note that if 



2N+1 

, fc+1/2 



P-W = - E (55) 



k=2 



then the formula (^) gives the realization of the general Krichever prescription for con- 



structing the algebraic-geometrical solutions to the Whitham equations |12 



4 Functionals of the Peierls-Frohlich type as the Poten- 
tials for the Local Solutions to the Whitham Equa- 
tions 

We point out the important relationship between the differential Aj and the standard mero- 
morphic differential (quasimomentum) 

dp = d\ (56) 



normalized by 



As can be easily seen. 



dp = 0, m = l,...,N. (57) 



A. = (5^ 
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In particular, for = 0: 

A = 2drdp = p^^^ (59) 



It is well known (see, for instance, |TT|) that / dp is the generating function for the averaged 
Kruskal integrals Qk- 

P^/..^2v/A.f^. (60) 
Expanding (|56D near infinity and comparing with (^O]) we find that 

I 2N+1 

qi = Qo-7^ J2 ' (61) 



where 



1 f 

Qo = u = lim — / u{x, r)dx (62) 



Then it follows from (|58|) that 



provided djU ^ 0. Substituting the representation (|63D into the solution (|3T1) one arrives at 
the potential form of the local solution to the Whitham equations 

a,F;v(r;a;,t) = 0, j = 1, . . . , 2Ar + 1 , (64) 

where the potential 



F;v(r;a;,t) = $[Mo(a^);a;,t,iV] = + 4tA^/2}rfp + 2 ^ p^{\)dp-2i J T{\)dp] , 

(65) 

In the case A (E=l) we observe that represents the functional of the Peierls-Frohlich 
type 0, Igip, 0: 



Fj^{r;x,t) = ^^[j^{xX^/^ + 4:tX^/^}dp + 2 J P-{X)dp)]. (66) 

Qoo UGfe 

The case of the Peierls-Frohlich type functional with (cf. ([55|)) 

P-(A) = - E c,X'^'^' + 2 -pl=du, (67) 



k=2 
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where g{u) is a sufficiently small smooth function, has been studied recently by Dubrovin 
[|19|1 who found that the minimizer to (H), (|67|) gives the solution to the Cauchy problem 
for the KdV- Whitham system with the monotone initial data 



2N+1 



^ = [ E + ^(«)]*=o • (68) 



{2k + 1) 

k=2 

One can suppose that the minimizer to the functional ( |65|) gives rise to the solution to the 
Cauchy problem with hump-like initial data. We emphasize also that the solution in the 
form (|64|), (|65|) does not require analyticity from the initial data. 



5 Global Solutions to the Whitham System 
5.1 General Formulation 

As it was shown in Sec. 2 which gives an account of the results of [|1|, the steepest descent 
for the RH problem yeilds the global solution to the Whitham equations bypassing the 
procedure of integration of the modulational system itself. The obtained solution, however, 
is restricted by the requirements of decaying at infinity and of analyticity imposing upon 
the initial data. In addition, the consideration in concerns only with solitonless (pure 
reflective) initial data. On the other hand, the Whitham equations themselves admit more 
general formulation of the problem cancelling these restrictions. 

Our idea is to construct the global solutions to the general Cauchy problem for the 
Whitham equations by applying the results of the RH problem approach described above 
to a more general class of functions p(A), t(A) (and therefore to a more general class of 
initial data) appearing in the solution (|3lD . More in detail, the functions p(A), r(A) can be 
multivalued (even infinite- valued) with different number of branches for different A. Actually, 
as we will show, such type of behavior for p and r corresponds to a multihump (or infinite- 
hump nondecaying) initial data. Also, the resulting formulation of the RH problem does 
not require any analyticity from the functions p(A) and r(A) which cancels requirement of 
analyticity for the initial data and is consistent with the hyperbolic nature of the zero- 



dispersion KdV limit |30 



First, we define following Dubrovin the Whitham system as a sequence of the mod- 
ulation systems (0) defined for = 1,2,... For A^ = this coincides with the Riemann 
wave equation 

dtr — Qrdxr = , (69) 

with the initial data r{x, 0) = uo{x) are given. Solutions of the Whitham equations (|34D for a 
given A^ typically exist only within certain domains of the [x, t) plane. The main problem of 
the theory of the Whitham equations is to glue together these solutions in order to produce 
a C^-smooth multivalued function of x that also depends C^-smoothly on the parameters 
ri, . . . ,r2Ar+i. 

Thus, to produce the global solution to the Whitham system with the aid of the obtained 
local solutions (^) one should: 
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i) determine the right genus at every point xo,to- 

ii) provide C^-smooth matching of the solutions (|3l|) for different genera on the phase 
transition boundaries [|lT], |]19 . 



The properties of the phase transitions can be investigated inside the local Whitham 
theory. 



5.2 Phase Transitions 



We study what happens to the solution (|3T|) when one of the Riemann invariants r2j coalesces 
either with r2j-i or with ^2^+1. It follows from the solution (PTf) that its phase transition 
properties are completely determined by the properties of the differential Aj near the double 
points Tfc = Tfc+i and can be investigated directly using the representation (^) for A^. It 
is more convenient, however, to use the relationship (^) between and the meromorphic 
differential dp properties of which are known well. 

Near the double points, the differential dp as well as the coefficients Qk in the decompo- 
sition (|60[) have the following asymptotics [19|: 



ij r2j+i-r2j 
fN{ru...,r2N+i) 

Here 



^ (small gap) 
fN-i{r I,. . . ,r2j,r 2j+i, 



,r2N+i)+J^^fh,i{ri,---,r2, 



'31 



faj+i,...;/?, i/)+o(z/2) 
(70) 



r2j+i - r2j 



(71) 



Here and below the hat means that the correspondent coordinate is omitted. 
2) r2j — T2j-i (small band) 



In in,. 
Here 



r2N+i) = /iv-i(ri, . . . , hj-i, hj, r27v+i)+5/^j(ri, . . . , hj, ■■■■,V, S)+o{S) . 

(72) 

V="-^'^^, 5= [log 7 ^-^r- (73) 



2 ' ^ (r2j -r2j_i)2 

One can see that in both cases, in the limit, the double points drop out of the function 
Jn and it turns into its analog for the N — 1 genus. This fact follows from the normalization 
for the meromorphic differential 

r2j+l 

I dp = 0, j = l,...,N, (74) 



N 



which implies that the polynomial \ + Qk^ in the numerator of dp has exactly one 

k=l 

zero in each gap. Then, if one shrinks either gap or band, this zero inevitably coincides with 
the double point in the denominator. 

Certainly, all the averaged Kruskal integrals Qk{ri, . . . ,r2N+i) (see (|60D) have the same 
properties near the double points. Then using the representation (^) for the generating 
differential we arrive at the following asymptotics by differentiating ( [7D| ) and (|72p. 
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and 



1) r2j+i - r2j 

^rVi, • • • , ^27v+i; A) = At^~^'(ri, ...,r2j, hj+i, • • • , r2iv+i; A) + 0(z/2) , (75) 

if k^2j,2j + l, 

APin,...,r2N+i;\) = A[';i + 0iu), (76) 

if = 2j or = 2j + 1 , 

where aI^^ is the hmiting value of the differential A^^^ which follows from (^), ( PBf ) when 
one pinches the j-th gap : 

_|_ 

M^l = ^f\n, r2j-i, r2j, r2j,r2j+2, • • • , raTv+i; A) = "' d\ , (77) 

{^-r2jY\ R2N-1 



R2N-1 = (A - n)(A - ra) ... (A - r2j_i)(A - r2j+2) ... (A - raTv+i), 

and Ag^l has a double pole at A = r2j. 

2) r2j - r2j-i 

• • • ,'"2jv+i; A) = Ajf^^Vi. • • • ,^2j-i,r2j, . . . ,r2Ar+i; A) + 0(5) , (78) 
if A;^2j-l,2j, 

and 

Afk^i, • • • ,^2Ar+i; A) = Aif + 0(z./5) , (79) 
if /c = 2j — 1 or k = 2j , 
where the Ag^l is the limiting value of the differential A\^^ when one pinches the j-th band : 

\N-l _|_ 

^2j- = ^^k^\fi, . . . , r2j_2, r-2j, r2j, rs^+i, . . . , r27v+i; A) = ^ dX , (80) 



(A-rajOy'i^^jv-i 
^2Af-i = (A - n)(A - r2) ... (A - r2j-2)(A - r2j+i) ... (A - r2Ar+i), 

and A^Jl has a single pole at A = r2j (cancellation of one pole occurs due to the zero in the 
vanishing j-th band ). 

One can see the substantial difference in the limiting behavior of the differential Ag^^ 
depending whether one shrinks the j-th gap ([77| ) or the j-th band (|80|). 

The asymptotics ( [75D and ( [75| ) provide natural - smooth matching of the solution 
( |3TD for different genera on the phase transition boundaries. The limiting values in (|76D 
and (^) determine the motion of those boundaries. Namely, the boundaries Xj± separating 
A^-phase and — 1-phase regions correspond to closing either the j-th gap (+) (linear wave 
degeneration) or the j-th band (— ) (soliton degeneration) and obey the ODE's : 



^ = -4i?e...A^/^A[;;^i 
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^1) 



One should also check what happens to the solutions (|3l|) under the phase transition 
(N = 1) ^ (N = 0). As = we have Af = (A - r)~^/^dX (see (H)) and, as a resuh, 
we arrive at two different equations corresponding to the cases A and B respectively. In the 
case A we have 

2t / 2d} p(X) „ , , 

floo ^ ^ 



which is the solution of the Riemann wave equation ( p9D with the initial data in the form of 
the increasing part of the hump Uo{x): 

X + 6rt = x~{r) . (83) 

For the case B we have an analogous result for the decreasing part of the initial hump uq : 

2d} t(X) 

B: x + 6rt = x'(r) / / ^ c^A = x+(r) . (84) 

r 

It can be easily seen also that the decompositions of A^^^ near A^°^ have the form ([TBI), 
( ffSD that provides - smoothness of the transition A^ = 1 ^ A^ = 0. This transition was 



investigated in detail by Avilov, Krichever and Novikov |^ (see also |TT|, [|T^]). 

It is clear now that the cases A and B in the solution ( pT]) describe different parts of the 
global solution of the Whitham equations corresponding to the monotonic branches of the 
initial profile uo{x). We also note that formulas (|82D, ( ^4]) give the solution to the Riemann 
equation in terms of the KdV scattering data and were obtained in this form by Geogjaev 
in 131 



5.3 Inequalities: Determination of Genus 

As it has been said before, the necessary part of constructing the solution to the Cauchy 
problem is determination of the genus of the Riemann surface at every point x,t. In fact, 
this problem cannot be resolved locally in the frame of the obtained solution (|3T|) for the 
Whitham system. This point of view has not been established clearly in the literature so we 
will describe it more in detail. 

The genus (the number of oscillating phases) A^ is the "external" characteristics with 
respect to the local solutions to the Whitham system. To illustrate that more clearly, we 
present a simple example of ambiguity in the local determination of genus. 

In the Figure 1, there is a typical solution of one-phase averaged Whitham equations for 
the initial data with the only break point (for example, the cubic breaking in the Gurevich - 
Pitaevskii problem 0) depicted by a solid line in the region (a;_, a;+) while the three- valued 
solution to the Riemann wave equation ( |69l) is depicted by a broken line in the region (xi, X2), 



Xi < x_ < X2 < x+. Outside of those intervals both solutions coincide. 

One can see that in the vicinity of any point of an open interval (x2,a;+) both cases 
(7 = and g = 1 can be applied without any violations of the existence for the local solution. 
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The missing part of the information which enables to distinguish the unique solution in this 
interval is the break time tcrit (if t > tcrit then = 1 otherwise = 0) but this information 
is relevant to the global properties of the solution for case = 1. In other words, generically 
to make the right decision about the genus at t = to one should follow the solution all time 
until to and change genus after passing the critical points . The rigorous global predictions 
of the genus have been made so far only in a few simple cases. 

In the case of analytic initial data with the only break point there are results by Tian 
[14], |15| that the solution beyond the break time globally belongs to = 1. 

Another result in this area is due to Grava and it reads that the maximal genus in the 
Cauchy problem with polynomial initial data does not exeed the degree of the polynomial 
and N —>■ 1 asymptotically as t — > cxd (the latter result also can be found in p). 

But even in the simplest case of the polynomial data the determination of genus at a 
fixed point Xo,to requires knowledge of the behaviour of the solution in all earlier times. 

The RH problem approach makes it possible to determine the right genus of the problem 
locally with the aid of the inequalities (H) - (|12D for the complex phase function g{X). These 
original inequalities were derived for the decaying solitonless initial data in It is 

clear, however, that if one assumes the finite speed of propagation (hyperbolicity) in the zero 
dispersion limit (which has been rigorously proved for the case of decaying initial potential 
0) [0! [0) then one can consider an arbitrary initial profile and regard these inequalities 
as a complementary part to the general local solution ( PTf ) for finite t. Really, given the 
functions p±(A) and r(A) as the simple integral (Abel) transforms (|3|),(|47|), (^) of arbitrary 
initial data one can construct the function g{X) and its derivative g'{X) and then, by trial 
and error, to determine the genus with the aid of inequalities. Also, one should add the 
relations @ and ( p!0D distinguishing the cases A and B for the problem with nonmonotone 
initial data. Certainly, for the monotonically increasing data we always are in the case A. 

Example: One-Hump Problem 

Following we present now an example of a systematic procedure for obtaining A^ = 
A^(a;, t) for all x at any given t. The result of this procedure is to construct for each t the 
separatrix F{t) = {{x,r) : r = r{x,t), —oo < x < oo,r = ri, . . . ,r2N+i} for the case of 
one-hump initial data with finite number of break points (see Figures 2a and 2b). 

For times t less than some critical value tcrit one takes A^ = A^(x, t) = and solves (|69D 
for r; the time tcrit corresponds precisely to the time at which (|I]) with e = breaks down. It 
turns out that the associated function g constructed as above indeed satisfies the auxilliary 
inequalities (|^), (||), (^, (|T0|) and hence g is the desired solution. Thus, for t < tcrit, we obtain 
a separatrix of the shape of Figure 2a; moreover, we find that F{t) = {x,r) : r = u{x,t), 
where u{x,t) is the solution of (|ip with e = 0. For x > Xo(t), we are in case B and for 
X < Xo{t), we are in case A. When t > tcrit, one again sets A^ = N{x,t) = for x >> 1, 
computes r from ( pPj ) and verifies once again that the side conditions are satisfied in case 
B. However, as we move x to the left, we find that at least one of the inequalities (^, (§), 
(H), ([ToD breaks down. In Figure 2b below the first inequality in (H), — r < ((?_!_ — g^)/2i, 
breaks down a.t x = Xi and for x < Xi the interval Gq breaks up into two intervals Go and 
Gi- For such x one solves the system ( pT]) for one of the cases A or B to get ri, r2, and 
and verifies that indeed the associated A^ satisfies the above inequalities. In the scenario of 
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Figure 2b, as we move x further to the left, the same inequahty for r again breaks down at 
some X = X2{t), for some r G Go{x2{t),t). Again Gq sphts up into two intervals yielding a 
total of three intervals Gq, Gi and G2, etc. As we continue to move towards X3, the interval 
Gi shrinks to a point and eventually disappears. For x between X4 and X3 we are again in 
the two interval case and for x < X4 we return to the one interval case N = 0. For x > Xq, 
we are always in case B and for x < Xq we are always in case A. The scenario of figure 2b 
is representative of the generic case in which any finite number of intervals may be open at 
some X. For more complicated, indeed pathological initial data infinitely many folds may 
appear. 

6 Zero-Dispersion KdV with General Initial Data 

The basic assumption which is made in the following consideration is the hyperbolicity of 
the zero dispersion limit for the KdV with an arbitrary initial potential. By applying this 
assumption one can extend the formulation of the RH problem made in Sec. 2 for the case 
of the solitonless alanytic initial perturbation to the case of an arbitrary initial data and to 
construct the solution to the corresponding Cauchy problem for the KdV-Whitham system. 

It is clear that due to the (assumed) hyperbolicity of the problem only a finite number 
of 'humps' in the initial perturbation will be involved into nonlinear interaction at any 
finite time at any particular point. Therefore, the solution to the problem for the Whitham 
equations with an arbitrary (nondecaying) initial data can be obtained by solving the "multi- 
hump" (decaying) problem in the zero dispersion limit and then by tending the number of 
humps to infinity. Certainly, the solution obtained in this way will be valid during limited 
time interval (presumably for t <^l/e) until the semiclassical spectrum (Riemann invariants 
of the Whitham equations) begins to compete with the fine spectrum of the initial potential 
and the problem loses its hyperbolic character. 

The main technical obstacle to the direct extension of the formulation (^ - ([121) to the 
case of a multi-hump initial data is the multivaluedness of the functions p+(A) and r(A) 
®y (i) appearing in the RH problem. To manage this difficulty it is convenient to first 
reformulate the RH problem for one-hump initial data introducing the finite reference point 
instead of the reference point at infinity (scattering from the right). 

6.1 Reformulation of the One-Hump Problem for the Finite Ref- 
erence Point 

We start with the function 



which, as can be checked directly, solves the RH problem (^ - ([T2|) for t = 0, = 0. 



00 




(85) 



X 
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We introduce a new function 

XO 

H{X, X, xo) = J (A^/^ - (A - uoix')Y/^)dx' , (86) 

X 

where Xq is an arbitrary fixed reference point. Without loss of generahty we put Xq at the 
maximum of Uo{x) (see Fig. 2a). 
One can see that 

g±{X,x,0) = H±- X^/^xo + p+Tir+, (87) 

where 



(A) = / {uo{x) - Xy/^dx . 



Then we have the relationships (at t = 0) 



^ = ^ X' XO + P+, (89) 



9+ -9- H+-H^ 
= r4 

We also introduce 



(90) 



r_(A) = r+(A) - r(A) = / {uo{x) - Xfl'dx . (91) 



Then the function H{X, x, xq) is the solution of the RH problem which follows directly from 

m - (0) ■■ 

In the region x < x (see Fig. 2a) which lies on the left from the hump and contributes 
to the spectral band Bi the correspondent relations follow from (^ (case A) taking into 
account (|8|), (|90|): 

T--- ^<-"- '^^) 

where 

ao{X,x) = X^^'^ix ~ Xo) . (93) 

Analogously, for the region x' < x < x^ lying under the hump and contributing to the gap 
Go we have from (^, (^: 

\ - =-^,. ^~<^, <^+- (94) 

And finally, in the region x > x'^ which lies on the right from the hump and again 
contributes to the band Bi (but it is the case B now) in the original RH problem (see (|T0|) ) 
we have 

Ha- — H H, -\- H , . . 
= ^+' 2 (95) 
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Also we have outside the interval (0, 1) from ([TT|), (|T2[) 



H+ + H_ = if A < 0, (96) 

H+-H_=0 if A > 1 . (97) 

Due to the finite speed of propagation in the zero dispersion KdV with decaying initial 
data the topology of the picture (Fig. 2a) does not change under the t-evolution (see Fig 2b). 
Finitely many folds which appear at some finite t determine the band-gap structure at any 
point X so that the region x < a;" on the {u, x) - plane of the initial data will provide the 
bands Bn, {n = 1, . . . , A^ + 1) (-Ba^+i corresponds here to the case A), the region x~ < x < x~^ 
will provide the gaps Gk, {k = 0, . . . , N), and the region x > x~^ always will correspond to 
the band -Bat+i in the case B. Then one can formulate the time-dependent RH problem for 
the finite t by the substitution reflecting the linear temporal evolution of the phase a(A) (|): 

ao(A, x) ao(A, x, t) = X^/^{x - Xq) + AX^^H (98) 



in the formulas (|9^) - 

Thus, we have effectively removed the infinite reference point from the formulation of the 
RH problem corresponding to the evolution of the one-hump perturbation. 

6.2 Evolution of an Arbitrary Profile 

The strategy of investigation of the evolution for an arbitrary initial perturbation uo{x) 
is essentially the same as in the previous subsection: we formulate the RH problem for 
the function H{X, x, Xq) (^) at t = starting with the fixed reference point Xq and then by 
assuming the finite speed of propagation (hyperbolicity of the zero-dispersion KdV) arrive at 
the time-dependent if-function by the substitution ( PE| ) for the phase ao ■ As was mentioned, 
the situation, however, is complicated by the multivaluedness of the functions p(A) and r(A) 
appearing in such a problem. 

We begin with the multihump initial potential and split up the strip < m < 1 in the 
u,x - plane of the initial data into two type of domains (see Fig.Sa): 

r - domains correspond to humps in the intitial perturbation and contribute to the gaps 
in the semiclassical spectrum. 

p - domains correspond to wells in the initial perturbation and contribute to the bands 
in the semiclassical spetrum. 

We fix the reference point Xq at some maximum of the function Uq{x) and enumerate the 
monotonic parts of Uo{x) for x > Xq labelling them by the roots of the equation A = Uo{x) 
(to avoid unnecessary notation complexities we put the number of humps on the right from 
the point Xq equal to the number of humps on the left): 

X-M ^ -^-M+l < • • • < < 1^ ^ < X-^ < . . . < X]^j_i < X~^.^ , (99) 

MeN, ^>0, ^<0. 
dX dX 
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We denote the domains of definition for eacli x^{X) as D^. Then for each A = A* one can 
define two subsequences: {a;~ .(A*)} and {x^.{\*)}, {mj = ±1,±2, ... ,±Mj , j = j(A) = 
1, . . . , 2M — 1 , x^. > such that only those participate in {x^.} for which A* G -D^ 

(FigS.b). 

For each j = j(A) one can then define the sequences {pm^}, {Tmj}'- 

rrij rij 

Pm,{X) = Re J {X-uo{x')y/^dx', Tn^{X)=Im j {\ - u^{x')Y'Hx' (100) 

Then the function H{X, x, xo, 0) defined by ( ^6l) and corresponding now to a multi-hump 
analytic initial data uo{x) is related to the original function g{X,x, 0) by (cf. (|87|)) 



if < X < , mj>0: g± = H± - X^^'^xq + p+ - pm,-i T i^m^ , (101) 

where p+(A) is redefined for the multi-hump perturbation in the following way: 

oo 

p+(A) = <A^/2 + / [AV2 - (A - uo{x')y/']dx' , (102) 



Here one should put po = and Xq = xq. If mj < then the relationship ( |101| ) is valid 
within the interval x^. < x < x^^_^_i. 

Then, at t = we have instead of (pOf), ( pO] ) within the mentioned intervals covering the 
entire axis: 



2 2 - -V . r-r 

9+ -9- H 



X'/^Xo + p+-pm,-i, (103) 

(104) 

Now we can formulate the RH problem which is satisfied by the function H. We substitute 
( |103| ), ( p.04| ) into the original RH problem (|^) - ( p!2[ ) taking into account new definition of 
P+(A) O: 

In each r - domain (x^_i < x < x^. for rrij > and x^. < x < x^^_^,i for < ) one 
has from (|]), (|]) 

H'+H'_,, H+-H_ , , 

= -a^ + p„^._i , r„^,_i < < Trr,^ . (105) 

In each p - domain < x < x^. we have from (|^, (||), (|T0|) 

= r^^ , -ao + P™,-i < ^ < -tto + Pm, • (106) 

We note that inequalities (||), (P), ( |10|) have converted into one double inequality. One should 
also put To = 0. 
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Outside the interval (0, 1) one has analogously to (p6D, (p^): 

H+ + H^ = if A < 0, 



(107) 



H+-H^ = if A > 1 . (108) 

As well as in the previous subsection, the time-dependent RH problem is obtained from 
( [I05D - ( iroq ) by adding the term iX^^H to ao{\, x) (see (0)). However, to retain validity of 
the formulated RH problem after this substitution one should assume hyperbolicity of the 
zero-dispersion limit for the KdV with general initial data. In this case the toplogy of the 
separatrix would not change under the evolution. 

Recalling that r - regions contribute to the gaps Gk , {k = 0, . . . , N) and p - regions con- 
tribute to the bands Bn ,n = 1, . . . , A^+1, one immediately obtains the function if (A, x, xq, t) 
and its derivative (cf. (1T9|) , ([l5|) ). For the sake of simplicity we write down the result as- 
suming that the initial data are such that each p - domain contributes to only one band 
(analogously, each r - domain contributes to only one gap). Consideration of more general 
data would only complicate notations adding nothing to substance. 

(109) 

. . rz 77^ ^( f '2p*M - 2^0 (/i, X, t) - VL^ dp f -2iT^{p) dp 

ii(A,x,Xo,tj = yit2Af+i(Aj 2^ / / , ^ 77—+ / 



Vi V ^2^+1 (/^ - ^) '''' bU V - ^) ' 

(110) 

where 



= -2{x - xo) / A^/V, - 8t / A^/V, +4j2ilp*m^j-^ I , (m) 

O-oo 0,00 ^ ^ Gm Bm-\-l 

fio = 0, p* = 

. Here 

{p:;(A)} C {pu,{\)] , {r;(A)} C {r,^.(A)} , m = 0, . . . , iV . (112) 

The algebraic equations determining dependence of the band-gap structure on x and t are 
obtained from the condition (^) djVLi = and can be eventually represented in a potential 
form (|6^ ) with the generalized Peierls - Frohlich type functional as a potential (see Sec. 3, 4) 

djF^{r;x,t) = 0, j = 1, . . . , 2iV + 1 , (113) 



M^;^,t) = ^^[f{{x-Xo)X'/' + 4tX^^'}dp-2Y.{Jpldp-z I r^dp)]. (114) 

loo ^ Gm Bm + 1 
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The system ( |113|) , ( |114|) represents the general local solution to the Whitham equations. 



Due to the properties of the meromorphic differential dp discussed in Sec. 5. 2 this solution 
provides - smooth matching on the phase transition boundaries where genus N of the hy- 
perelliptic surface changes. The right choice of the local genus and the subsets {r^}, {Pm} 
( |112| ) is verifed by check of the fulfilment of the inequalities ( |105| ) and ( |106| ) for the function 
H{\, xq, X, t). 

As an example of an effective choice of the subsets {t^}, {Pm} '^^^ can indicate the 
case of periodic initial data uo{x) where = {m + l/2)/m J^^^'^iX — uo{x')Y^'^dx' , p*^ = 
mRe /^Q°^^(A — uo{x')y^'^dx' , T is the period, m G N 

Thus, we have bypassed the inverse scattering transform for the zero - dispersion KdV 
with general initial data by assuming hyperbolicity of the zero dispersion limit. One can see 
that the resulting construction requires neither decaying nor analyticity for the initial data. 
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Appendix 1: Proof of the Consistency Lemma 

Consider the function 

(J) _ djQj 



We rewrite the expression ( ]20|) in the form containing contour integrals. For analytic 
functions p+(A) and r(A) we have 

= {-2xX^^^ - StX^/^jipj + 2p+{X)tfjj - 2i j r{X)il)j. (Al.2) 

Here we have used the fact that the functions p+(A) (^) and r(A) @) have their own different 
Riemann surfaces with the branch points: 

0, oo for P+(A) 

1, oo for r(A) 

Then we obtain from (|A1.1|) 



$S^-) = -2/{xA^/2 + 4tAn|t' + 2 /p+(A)||^-2^ / r{xf^ . (Al.3) 
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We consider the differential 



Substituting (^21) into ( |A1.4| ) we immediately obtain 

(A - ri)yJR2N+l 

We integrate K^P over the a-cycles. Then, taking into account (^31), (|A1.4|) we get 

A^^ = 0, (A1.6) 



m,j = l,...,iV, z = l,...,2iV + l 
The normalization (|A1.6|) uniquely defines the coefficients pk,i independently on j, i.e. 

A«=Af) = ... = Af)^A,. (A1.7) 

Therefore, the differential ( [A1.4| ) does not depend on the index j, and, according to ( |A1.3| ), 
the function does not depend on j as well. That means that the system (|30D is cosistent 
and takes the form 

j{x\^'^ + At\^/^}kj= j p+{\)kj-i j r(A)Aj-, (A1.8) 



provided dikj ^ 0. 
Q.E.D. 



Appendix 2: Derivation of the Identity (^) 

We introduce the analytic function 

oo 

/(A) = j [Ai/2 - (A - u,{x')f/^]dx' = p_(A) + p+(A) - 2r(A) . (A2.1) 

— oo 

Then one can observe that the following identity holds 

ffiX)A, = 0. (A2.2) 

Really, /(A) ~ A"^/^ and A^- ~ X'^^'^dX as A — oo. Considering ( [A2.2| ) as the integral over 
the contour surrounding the interval (0, 1) we arrive at the identity for analytic P-(A), p+(A) 
and r(A). 

^ {p_(A) + p+(A)}A, - i j r(A)A, = . (A2.3) 



24 



References 

[1] P.Deift, S.Venakides, and X.Zhou, New results in small dispersion KdV by an exten- 
sion of the steepest descent method for Riemann-Hilbert problems, International Math. 
Research. Journ., No 4 (1997), 285-299. 

[2] G.A.El, Generating function for the Whitham-KdV hierarchy and effective solution of 
the Cauchy problem, Phys.Lett.A 222 (1996), 393-399. 

[3] A.V.Gurcvich and L.P.Pitaevsky, Nonstationary structure of a collisionless shock wave, 
Sov.Phys.JETP 38 (1974), 291-297. 

[4] G.B.Whitham, Non-linear dispersive waves, Proc. Roy. Soc. A283 (1965), 238-291. 

[5] H. Flaschka, G. Forest, D.W. McLaughhn, Multiphase averaging and the inverse spectral 
solutions of the Korteweg - de Vries equation. Comm. Pure Appl. Math. 33 (1979), 739- 
784 . 

[6] P.D.Lax and C.D.Levermore, The small dispersion limit of the Korteweg - de Vries 
equation I, II, III, Comm. Pure Appl. Math. 36 (1983) 253-290, 571-593, 809-829. 

[7] S.Venakides, The zero dispersion limit of the Korteweg - de Vries equation for initial 
potentials with nontrivial reflection coefficient. Comm. Pure AppL Math. 38 (1985), 
125-155 . 

[8] S. Venakides, The generation of modulated wavetrains in the solutionof the Korteweg - 
de Vries equation. Comm. Pure Appl. Math. 38 (1985), 883-909. 

[9] S. Venakides, The zero dispersion limit of the Korteweg - de Vries equation with periodic 
initial data, AMS Transaction 301 (1987), 189-226. 

[10] S.P.Tsarev, On Poisson brackets and one- dimensional systems of hydrodynamic type, 
Soviet Math. Dokl. 31 (1985), 488-491 . 

[11] B.A.Dubrovin and S.P.Novikov, Hydrodynamics of weakly deformed soliton lattices. Dif- 
ferential geometry and Hamiltonian theory, Russian Math. Surveys 44 (1989), 35-124 

[12] I.M. Krichever, The method of averaging for two-dimensional "integrable" equations, 
Funct. Anal. Appl. 22 (1988), 200-213. 

[13] I.M. Krichever, The T-function of the universal Whitham hierarchy, matrix models and 
topological field theories, Comm.Pure Appl.Math., 47 (1994), 437-475. 

[14] F.R.Tian, Oscillations of the zero dispersion limit of the Korteweg - de Vries equation. 
Comm. Pure Appl. Math. 46 (1993), 1093-1129. 



25 



[15] F.R.Tian, The initial value problem for the Whitham modulation equations, Comm. 
Math. Phys. 166 (1994), 79-115 . 

[16] A.V.Gurevich, A.L.Krylov and G.A.El, Riemann wave breaking in dispersive hydrody- 
namics, JETP Lett. 54 (1991), 102-107 , Evolution of a Riemann wave in dispersive 
hydrodynamics, Sov. Phys. JETP 74 (1992), 957-962. 

[17] A.V.Gurevich, A.L.Krylov, N.G.Mazur and G.A.El, Evolution of a localized perturba- 
tion in Korteweg - de Vries hydrodynamics, Sov. Phys. Dokl. 37 (1992), 198-200 , 
A.L.Krylov, V.V.Khodorovskii, G.A.El, Evolution of a nonmonotonic perturbation in 
Korteweg - de Vries hydrodynamics, JETP Lett. 56 (1992), 323-327 . 

[18] T.Grava, On the Cauchy Problem for the Whitham Equations, Ph.D Thesis, SISSA 
(1998). 

[19] B.Dubrovin, Functionals of the Peierls - Frdhlich type and the variational principle for 
the Whitham equations, Amer.Math.Soc.Transl.(2) 179 (1997), 35-44. 

[20] E.D.Belokolos, Peierls - Frdhlich problem and potentials with finite number of gaps I, 
II, Theor.Math.Phys, 45 (1980), 1022-1026, 48 (1081), 604-610. 

[21] I.M. Krichever, The Peierls model, Func. Anal. Appl. 16 (1982), 248-263. 

[22] S.A. Brazovskii, I.E.Dzyaloshinskii, I.M. Krichever, Discrete Peierls - Frdhlich models 
with exact solutions, Sov.Phys.JETP 56 (1982), 212-225. 

[23] A.B.Shabat, One dimensional perturbations of a differential operator and the inverse 
scattering problem. Problems in Mechanics and Mathematical Physics, Nauka, Moscow, 
1976. 

[24] B.Beals, P.Deift, and C.Tomei, Scattering and inverse scattering for first order systems, 
AMS, Math.Surveys and Monographs No.28 (1988). 

[25] S.Venakides, The small dispersion KdV equation with decaying initial data, Proc.Symp. 
Appl. Math 54 (1998), 369-392. 

[26] V.R.Kudashev, "Waves-number conservation" and succession of symmetries during a 
Whitham averaging, JETP Lett. 54 (1991), 175-178. 

[27] P.D.Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of 
Shock Waves, Conference Board of the Mathematical Sciences, Regional Conference 
Series in Apphed Mathematics, SIAM, Philadelphia (1973). 

[28] G.V.Potemin, Algebro- geometric construction of self-similar solutions of the Whitham 
equations, Russian Math. Surveys, 43 (1988), 252-253. 

[29] O.Wright, Korteweg - de Vries zero- dispersion limit: through first breaking for cubic-like 
analytic initial data. Comm. Pure Appl.Math., 46 (1993), 421-438. 



26 



[30] C.D.Levermore, The hyperbolic nature of the zero dispersion KdV limit, Comm. Partial 
Differential Equations 13 (1988), 495-514. 



[31] V.V.Avilov, S.P.Novikov, Evolution of Whitham zone in KdV theory, Sov.Phys.Dokl. 
32 (1987), 366-368; V.V.Avilov, l.M. Krichever, S.P.Novikov, Evolution of the Whitham 
zone in the Korteweg - de Vries theory, 564-566. 

[32] V. V. Geogjaev, On the quasiclassical approach in the inverse scattering transform 
method, Dokl. Akad. Nauk SSSR 284 (1985), 1093-1097 (in Russian). 

Figure Captions 

• 1. Ambiguity in local determination of genus. Solid line: solution of the Whitham 
equations {N — 1), broken line: solution of the Hopf equation {N — 0). 

• 2. Separatrix {F{t) — {{x, r) : r{x, t) — ri, . . . , r2Ar+i}} , one-hump case, 
a) F(0) = {{x, r) : r = Uo{x)} b) F{t), t > t^Ht 

• 3. Multi-hump initial data, a) p- and r- domains, b) Sequence {x^.}. 
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Fig. 3a 
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Fig. 3b 



